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Abstract. In the past years modern mathematical methods for image analysis have 
led to a revolution in many fields, from computer vision to scientific imaging. How¬ 
ever, some recently developed image processing techniques successfully exploited by 
other sectors have been rarely, if ever, experimented on astronomical observations. 
We present here tests of two classes of variational image enhancement techniques: 
’’structure-texture decomposition” and ’’super-resolution” showing that they are effec¬ 
tive in improving the quality of observations. Structure-texture decomposition allows to 
recover faint sources previously hidden by the background noise, effectively increasing 
the depth of available observations. Super-resolution yields an higher-resolution and a 
better sampled image out of a set of low resolution frames, thus mitigating problemat¬ 
ics in data analysis arising from the difference in resolution/sampling between different 
instruments, as in the case of EUCLID VIS and NIR imagers. 


1. Structure-texture image decomposition 


A general approach to the denoising problem is based on the assumption that an image 
/ can be regarded as composed of a structural part w, (i.e. the objects in the image), and 
a textural part v which corresponds to finest details plus the noise. In this paper we have 
considered image decomposition mode ls based on t otal v ariation regularization meth¬ 
ods following the approach described in lAujol et al.l (l2006h . Such image decomposition 
technique is based on the minimization of a functional with two terms, one based on the 
total variation and a second one on a different norm adapted to the texture component. 
Given an image / defined in a set Q, and being 5y(Q) the space of functions with lim¬ 
ited total variation in Q (i.e. V{f, Q) = f |V/(x)| dx < +cx)) we can decompose / into 


its two componentes by minimizing: 
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where ||.||^ indicates the norm of a given space X (e.g. L^(Q)) and the minimum is found 
among all functions (u,v) e (BV(Q) x X) such that u + v = f. The parameter is a 
natural exponent, and A is the so-called splitting parameter. The best decomposition is 
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Figure 1. From left to right: simulated noise-free mosaic, simulated observed 
mosaic (CANDELS-GOODS depth), ’’Structure” image obtained thorugh TVL2 de¬ 
composition. 


found at the A for which the correlation between u and v reaches a minimum. We have 
designed a C++ code, named ATVD (Astro-Total Variation Denoiser) which implements 
three versions of the technique, based respectively on the TV-L^ (X = L^(Q) ), TV-L^ 
(X = (Q)) and TVG (X being a Banach space as defined in lAujol et all (l2006l) ) norms. 

Tests on simulated HST images. We have tested the performance of the decom¬ 
position technique method on a 15x15 arcmin simulated image reproducin g the average 
depth of the F160W (H160 hereafter) CANDELS GOODS-South mosaic (iGrogin et al. 
I2OIII) . The analysis of a simulated image in which the position of sources below and 
above then nominal detection limit are known, allows to test the effectiveness of de¬ 
composition techniques at increasing the depth of astronomical observations. We de¬ 
composed the simulated image through our ATVD code using different norms and we 
found that the TV-L^ one is the most effective. The best Structure image (H160^^ here¬ 
after) obtained is shown in FigH] together with the input noisy one (H160/^p) and, as 
a reference, the simulated image without noise. At a visual inspection the resulting 
Structure image appears deeper and richer in details than the input observed mosaic. To 
qua ntitatively assess the adv antages of the TV-L^ decomposition we have run SExtrac- 
tor (iBertin & Amoutsl[l996l) on both H160^^ and H160/^p. By varying SExtractor input 
parameters affecting source detection we have estimated completeness and purity of the 
extracted catalogues. A comparison of the best detection performances on H160^^ and 
H160/^p is shown in EiglJl Source extraction on the Structure component H160^^ yields 
to an higher purity of the catalogue (~0.5 mags deeper confidence threshold) at a sim¬ 
ilar completeness reached on the original H160/^p. Alternatively, by pushing detection 
to lower S/N, we can obtain a much higher completeness with acceptable contamina¬ 
tion levels (right panel in EiglJ]) with respect to H160/^p. As shown in Eigl2]magnitudes 
estimated on H160^^ are reliable also for those sources up to mag~30 which are unde¬ 
tectable in the original mosaic. 


2. Super-resolution 


A low-resolution, blurred frame can be considered as the result of the application of a 
number of operators (downsampling D, warping W, convolution F/), to an ideal hig h 
resolution image, plus the addition of a noise component (e) fe.g. llrani & Peledri99lh : 
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Figure 2. Top panels: magnitude counts on the original and on the denoised 
(Structure) image maximising purity (left) or completeness (right) on the latter. Bot¬ 
tom panels: on the left, comparison between magnitudes measured on the denoised 
image and input magnitudes (additional faint sources not found on the noisy mosaic 
are shown in green). Right: sources detected on the Structure image but not in the 
input noisy one as they appear in the noise-free simulated mosaic, in the noisy image, 
and in the Structure image. 


= DHWx^ + c = Kx^ + e. Following this image formation model, it can be 
shown that given N low resolution frames with pixel sampling z, we can recover an 
higher resolution well-sampled image with pixel-sampling zj ^[{N) by minimizing the 
following energy function: 

E{x^) = \i\\^L - + m^H)) ( 2 ) 


where the last regularization term, which is needed for the problem to be well-posed, 
is a given norm of the image gradient. The minimization of Eq. [2|is found iteratively 


starting from a first guess on xh, e.g. a bilinear or polinomial interpolation of xl. We 
have designed a FORTRAN90 code (SuperResolve) to implement variational super¬ 


resolution reconstr^ion_of astr onomical images based on and regularization 
following lUnger et~^ (1201 Ol) and Zomet & Peled (l2002h respectively. 

Tests on simulated EUCLID images. The Euclid satellite will observe >15000 
sq. deg. of the sky with a visible (VIS) and near-infrared imager (NIR) to constrain 
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Figure 3. Top: simulation of a VIS image, one NIR frame, super-resolved NIR 
TV-L^ and TV-L^ images obtained from 9 NIR frames. Bottom left: PSF-FWHM 
measured on: the VIS image, the SR-NIR at the VIS pixel-scale, and on a NIR 
coadded mosaic resampled to the VIS pixel-scale. Right: comparison between the 
S/N of sources in the SR-TVL2 mosaic and in a standard coadded NIR mosaic. 


cosmological parameters to an unprecedented accuracy dLaureijs & et al.ll201lh . SR 
reconstruction can be of particular interest to match the resolution (sampling) of VIS 
(pix-scale=0.1”, PSF-FWHM=0.2”) and NIR (pix-scale=0.3”, PSF-FWHM=0.3”) im¬ 
ages thus improving reliability of photometric analysis. Through a dedicated simu¬ 
lation code we have produced simulated VIS and NIR frames by degrading existing 
HST-CANDELS images. We applied SR-L2 and SR-LI reconstruction to a set of 9 
NIR single epoch frames to produce a super-resolved mosaic at the same pixel sam¬ 
pling (0.1”) of the VIS mosaic. As shown in Fig. [3] SR-mosaic show an higher level 
of details, and the PSF-FWHM measured on them is not as worsened with respect to 
the real value as in the case of pixel upsamling through interpolation. The SR-L2 tech¬ 
nique is particularly promising since fluxes of the sources are preserved and the final 
S/N reached is the same as expected in a standard coadded mosaic. 
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